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A novel simulation strategy is proposed to search for semiconductor quantum devices which are 
optimized with respect to required performances. Based on evolutionary programming, a tecnique 
implementing the paradigm of genetic algorithms to more complex data structures than strings 
of bits, the proposed algorithm is able to deal with quantum devices with preset non-trivial con- 
straints (transition energies, geometrical requirements, etc.). Therefore, our approach allows for 
automatic design, thus avoiding costly by-hand optimizations. We demonstrate the advantages of 
the proposed algorithm by a relevant and non-trivial application, the optimization of a second- 
harmonic-generation device working under resonance conditions. 
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In addition to the new fundamental physics |-which 
can be investigated in lower-than-3 dimensions^, the 
enormous interest in semiconductor nanostructures has 
been strongly driven by the possibility of engineering 
semiconductor-based quantum devices (SBQD's)B; In 
particular, among the most successful applicatioas one 
should mention unipolar (i.e., intersubband) lasersO — the 
so-called quantum-cascade lasers — as well as a variety of 
semiconductor-based photodetectorgj. 

The design and optimization of SBQD's is often 
grounded on an expensive trial-and-error process: after 
a "promising" structure has been identified, an ad-hoc 
optimization is performed with respect to a few struc- 
tural parameters. In this Letter we propose a novel 
simulation strategy based on evolutionary programming 
(EP)l3, to efficiently search in a large parameter space for 
optimal SBQD's with respect to required performances: 
quantum-device optimization by evolutionary search (Q- 
DOES). 

The optimization of SBQD's is a complex multidimen- 
sional maximization problem of a fitness function con- 
taining the quality factors of the device; the hardness 
of the problem grows rapidly with the number of free 
parameters (number and width of the constituent lay- 
ers, alloy and dopant concentrations, etc.); furthermore, 
it may be required that the device works under certain 
conditions (e.g., with a high optical efficiency in a well- 
defined energy range) leading to a non-trivial constrained 
problem. In this respect, genetic algorithms (GAs)eI are 
often invoked for complex optimization tasks; in a typi- 
cal GA, a pool of possible solutions to the problem — the 
population — is evolved following a genetic paradigm: so- 
lutions belonging to the population — coded as strings of 
bits, called chromosomes — are allowed to mate in pairs 
in order to cross their genomes (genetic crossover), while 
random mutations (bit fiips) take place at a given rate 
to ensure genetic diversity. At the heart of the algorithm 
is selective pressure: the fittest chromosomes are given 
more chances to mate and pass their genomes to the next 



generation. 

For complex data structures the bit-string represen- 
tation is not usually the most natural one. In such 
cases, the genetic paradigm may be implemented more 
efficiently by designing context-dependent genetic opera- 
tors acting on non-conventional chromosomes represent- 
ing the "natural" data structure of the problem with im- 
proved efficency, a, s.trategy often referred to as evolution- 
ary programming^. In EP, genetic operators may have 
a direct physical interpretation, and may be optimally 
designed to exploit the physical insight of the problemQ. 
Finally, in EP inclusion of constraints is usually simpler. 

Let us consider a conventional multilayered het- 
erostructure formed by a fixed number, Ni, of layers of a 
PlxBi-x alloy; layer i is characterized by the A concentra- 
tion, Xi, and the number of A^Qi-x monolayers, n^. We 
map the structure onto a chromosome, each layer being 
in a one-to-one correspondence with a gene. Hence, we 
design each chromosome as a vector {gi, . . . , (?Ar, ), where 
gi contains all informations concerning a given layer; in 
the present case gi = (xi, ni)u. JDften, structural param- 
eters have physical limitationsLi; such constraints define 
the solution domain. It is important (and easy) to de- 
sign crossover and mutation operators which automati- 
cally produce offsprings within the solution domainllJ. 

We next discuss how to define the fitness function. For 
a SBQD the quality factors can be vpitten as a functional 
of the free-carrier wavefunctions V'n". From these ingre- 
dients, we can define a functional V[{'!/'n}] which mea- 
sures the performance of the device to be optimized; its 
evaluation may require the calculation of energy levels, 
optical matrix elements, scattering rates, etc., depending 
on the specific device requirements. Non-trivial (non- 
linear) constraints, such as multiple inter-subband tran- 
sition energies, may be dealt with by means of a penalty 
functional which is a measure over the parameter space, 
^[{■0"}] ^ 0, going to zero when all constraints are sat- 
isfied. In our implementation we always design V to be 
minimum for the fittest structures; more specifically, we 
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minimize 

V [{V„}] X exp {±V [{^„}] /5p{t)} ■ (1) 

5p{t) is a parameter which sets the strength of the con- 
straints and is, in general, a function of the simulation 
time t; +(— ) is taken in the exponential factor if V is 
positive (negative) definite, at least in the region of in- 
terest. 

The proposed Q-DOES algorithm proceeds as follows. 
We "evolve" a population of Np solutions, which is ran- 
domly initialized; at each simulation step we randomly 
choose Np < Np solutions which undergo, in pairs, dif- 
ferent types of crossover operations (see Tab. 1) with 
probabilities pc, designed so that if parents belong to 
the desired solution domain, so do the offsprings. Then, 
several mutation operators are applied (see Tab. 1), with 
probabilities to randomly chosen g^'s genes of the 
offsprings. Each mutation operator has a different pur- 
pose: e.g., boundary and uniform mutations ensure er- 
godicity, while non-uniform mutations perform a local 
search around the solution coded by the chromosomes 
which reach the later stages of the simulation. When Np 
offsprings have been generated by crossover and muta- 
tion, we convert the alloy concentrations and layer widths 
coded in their chromosomes — genotype — to the corre- 
sponding confinement potentials and space-dependent ef- 
fective masses — phenotype — and we evaluate the fitness 
function of the offsprings; in the present implementation 
the electron wavefunctions '0n and energy levels e„ are 
calculated with«3-|the conventional effective mass approx- 
imation (EMA)ll3. Finally, the set of Np + Wp chromo- 
somes are ranked according to their fitness-function val- 
ues, and the fittest Np chromosomes are passed to the 
next generation. 

Let us consider a technologically relevant and highly 
non-trivial optimization task: the design of a SBQD for 
second harmonic generation (SHG) working under reso- 
nance condition. Schematically, this is an equally spaced 
three-level structure, where the lowest inter-subband 
transition energies, £2 — ei — £3 — £2, are in resonance 
with a pumping radiation, huj, and light emission takes 
place at £3 — £1 = 2?itj_ At resonance, the emitted 
power is proportional tc&B /i — |/ii2/i23/^3i|, where fiij = 
J 'ip*{z)zipj(z)dz are the intersubband dipole matrix ele- 
ments corresponding to the individual optical transitions. 
Due to quantum confinement, these quantities can be 
greatly enhanced in a properly designed asymmetric het- 
erostructure with respect to bulk systems; however, to 
maximize p varying the confinement potential is a chal- 
lenging task, not only because the resonance condition 
should be preserved in the process, but also because /i 
is unbounded, and may become very large if resonance 
conditions are relaxed. Therefore, this example is ideally 
suited to test the proposed simulation scheme applied to 
a "hard" problem. This SHG optimization task has been 
undertaken so far only for extremely simplified-430tential 
profilegl3 and by means of analytical methodgi^. 



In order to maximize the SHG emitted power, we set 
V = —/X and we evolve a large population of type-I 
{x < 0.4) Al2:Gai_2;As heterostructures according to the 
Q-DOES algorithm described above. Resonance condi- 
tions with the CO2 laser radiation was enforced by set- 
ting V = [(£2 - £1 - hujf + (£3 - £2 - nujf]^/^. The op- 
timized structure produced by the algorithm is shown 
in Fig. 1: it corresponds to an estimated SHG power 
/i = 4.10 nm^, while the intersubband transition ener- 
gies fit hu) within numerical accuracy; this compares well 
with (and is actually 'v-f % better than) the largest esti- 
mated value 3.91 nm'^OJ obtained so far, to our knowl- 
edge, for the saHie hto, for an idealized continuously 
graded potentiaO. It is important to stress that the 
remarkably simple structure shown in Fig. 1 has been 
obtained with a random initial population, i.e., no guess 
has been made on the shape of the potential. 

The device in Fig. 1 is the best performing one in a 
small set of optimized structures corresponding to dif- 
ferent local minima the alghoritm converged to during a 
large set of runs, with corresponding fitness function val- 
ues falling within ~ 5% from the best value. The rate of 
convergence of our algorithm in a typical run is demon- 
strated in Fig. 2, where we show the values of the fitness 
function, Eq. (1), and of the SHG intensity, fi, of the full 
population. During the initial stages, the SHG intensity 
oscillates strongly, exploring a large part of the parame- 
ter space, until the search focuses on regions where the 
energy constraints are approximately satisfied. There, 
the exponential factor in (1) is ~ 1, and the SHG inten- 
sity coincides with the fitness function (in the figure the 
two vertical axes are shifted for clarity). Eventually, the 
distribution of the fitness function values gets narrower, 
as the population becomes degenerate, i.e., all chromo- 
somes become similar; at this stage, we always find the 
constraints to be satisfied to a very high degree. 

Although demonstrated with a simple, yet realistic 
problem, the algorithm is independent on the detailed 
structure of the genes, and more complex SBQD's can 
be dealt with using basically the same scheme; natu- 
ral candidates include quantum cascade lasers, mid- and 
far-infrared photodctectors, cavity resonators, etc. Be- 
sides, unlike recently proposed analytical optimization 
methodslij, a more accurate (microscopic) description of 
electronic states may be implemented using, e.g., tight- 
binding or pseudo-potential Hamiltonians. Obviously, 
dealing with more complex structures and/or Hamiltoni- 
ans may increase the computational cost of the algorithm 
(which was not an issue in the present implementation), 
which is dominated by the evaluation of the fitness func- 
tion; however, it should be noted that, in addition to the 
rapid convergence demonstrated above, the evaluation of 
the fitness function for different chromosomes is an in- 
trinsically parallel task which may take advantage of a 
parallel architecture. 
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typically require < Xi < 0.4. 

In the proposed implementation, we have inspired ourselves 
to the GENOCOP system described in Ref. |, but addi- 
tional operators have also been introduced which are suited 
to the specific problem 

One can write, at least in principle, any property of the 
quantum structure in terms of the orthonormal basis set 
obtained by the single-particle wavefunctions. 
Parameters are as follows: the Al^Gai-^^As band-gap Eg 
is obtained from E^ = Ef^^"" + 1.36x + 0.22x^ [C. Bo- 
sio, J. L. Staelhi, M. Guzzi, G. Burri, R. A. Logan, Phys. 
Rev. B 38, 3263 (1988)]; the valence band contribution 
is AE^ = 0.48a; [E. T. Yu, J. O. McCaldin, T. C. McGiU, 
Solid State Phys. 46, 2 (1992)]. The space-dependent efltec- 
tive mass is = 0.067-1-0. 083a: [Landolt-Bornstein, Semi- 
conductors: physics of goup IV elements, and III-V com- 
punds, O. Madelung, Ed. (Springer- Verlag, Berlin, 1982), 
vol. 17]. 
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The optimized potential in Ref. ^ spans more than 0.5 eV, 
and cannot be implemented as a strictly type-I structure. 
It is a distinct advantage of the present method the easy 
implementation of physical limitations on the parameter 
space, such as alloy concentrations. 



TABLE I. Description of crossover (C) and mutation (M) mechanisms 



C one-point 


parents' giS are split; offsprings t/i's are generated by recombination 
of the different parts 


C global arithmetical 


offsprings' gi's are obtained as linear combinations of the parents' (ji's 


M uniform 


gi is changed randomly with uniform distribution within the domain 


M non-uniform 


gi is changed randomly with a probability distribution which gets narrower 
as the simulation proceeds 


M boundary 


gi assumes values on the domain boundaries 


M average 


in gi and gj, Xi and xj are both set to {uiXi -f njXj)/(ni -\- nj) 


M flip 


gi and gj are interchanged 
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FIG. 1. Potential profile and for the lowest three confined states, shifted vertically by the confinement energy. The 

simulation was performed with Ni = 30 and Np = 100; Probabilities Pc and p,„ were in the range [0.07,0.1]. Resonance with a 
pumping radiation fiu) = 116 meV was enforced. The estimated = 4.10 nm^ corresponds to /ii2 = 1.725 nm, /i23 = 2.609 nm, 
/ttsi = 0.9104 nm. 
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FIG. 2. Fitness function values, Eq. (1) (top right axis), and SHG intensity, /x (bottom left axis), for a typical run (not the 
same leading to Fig. 1, but with same parameters). The two axis are vertically shifted for clarity. The solid line in the SHG 
intensity panel traces the fittest chromosome at each iteration. 
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